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VLF/LF  Reflectivity  of  the  Polar  Ionosphere 
20  July  — 20  September  1975 


1.  INTRODUCTION 

This  paper  provides  a summary  of  high  latitude  ionospheric  reflectivity,  as 

observed  by  the  USAF's  high  resolution  VLF/LF  ionosounder  operating  in  northern 
1 2 3 4 5 

Greenland.  * ' • * As  shown  in  Figure  1,  the  transmitter  is  located  at  Thule 
Air  Base,  Greenland  (76*  33'N.  Lat. . 68"  40'W.  Long.),  and  the  receiving  site  is 
106  km  north  at  the  Danish  Meteorological  Institute's  Ionospheric  Observatory  in 
Qanaq,  Greenland  (77*  24’N.  Lat.,  69'  20'W.  Long.,  Geomagnetic  Lat. 

89'  06'  N. ).  The  ionosounding  transmissions  consist  of  a series  of  extremely  short 
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Figure  1.  Geometry  of  the  Propagation  Path 


(less  than  150  peed  VLF  pulses,  precisely  controlled  in  time,  and  radiated  from 
a 130  meter  vertical  antenna.  At  the  receiving  site,  orthogonal  loop  antennas  are 
used  to  separate  the  two  polarization  components  of  the  ionospherically  reflected 
skywave  signal.  One  antenna,  oriented  in  the  plane  of  propagation,  is  used  to  sense 
the  groundwave  and  the  "parallel"  component  of  the  downcoming  skywave.  The 
second  loop,  nulled  on  the  groundwave,  senses  the  "perpendicular"  skywave  compo- 
nent. The  signal  from  each  of  the  antennas  is  digitally  averaged  to  improve  the 
signal-to-noise  ratio  of  the  individual  received  waveforms  before  they  are  recorded 
on  magnetic  tape.  An  example  of  the  observed  waveforms  is  given  in  Figure  2, 
where  the  parallel  waveform  (Figure  2a)conai8tsof(l)agroundwavepropagatedpulse,  i 

(2)  a quiet  Interval,  followed  by  (3)  the  first-hop  skywave  signal.  The  perpendicular 
waveform  is  shown  in  Figure  2b. 
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Figure  2.  Example  of  the  Observed  Waveforms 

Ionospheric  reflection  parameters  are  derived  by  computer  (AFGL's  CDC  6600) 
processing  of  the  ground  and  ionospherically  reflected  waveforms  with  allowance 
made  for  factors  such  as  ground  conductivity  (see  Section  4)  and  antenna  patterns. 

Although  the  data  are  recorded  about  once  per  minute,  for  this  paper  the  wave- 
forms are  averaged  into  2-hr  time  blocks  and  the  resulting  information  is  presented 
in  a weekly  format  (Figures  3 through  11)  as  described  below. 

2.  OBSERVED  WAVEFORMS 

In  part  A of  Figures  3 through  11a  set  of  averaged  parallel  and  perpendicular 
waveforms  is  presented  for  the  time  block  centered  near  local  noon  of  the  indicated 
day.  Each  of  these  waveforms  is  comprised  of  256  digitally  averaged  points  spaced 
2 Msec  apart.  In  part  B of  the  figures,  the  groundwave  Fourier  amplitudes  are 
shown  as  a function  of  frequency.  Although  the  data  presented  in  parts  C through  L 
of  the  figures  are  generally  limited  to  frequencies  in  the  first,  or  principal,  lobe  of 
the  spectrum,  information  at  higher  frequencies  can  be  used  when  sufficient  signal- 
to-noise  conditions  exist.  There  is,  however,  a frequency  range  around  each 
spectral  null  where  insufficient  signal  exists  for  measurements.  During  the  period 
covered  by  this  report,  the  transmitted  pulse  waveform  was  changed  to  lower  the 
frequency  of  the  maximum  spectral  amplitude  from  about  17  kHa  to  11  kHz. 
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3.  REFLECTION  HEIGHTS 


The  group  mirror  height  (GMH)  of  reflection  was  obtained  by  determining  the 
group  delay  of  the  skywave  relative  to  the  groundwave  and  attributing  this  time 

difference,  by  simple  geometry  (assuming  a sharply  bounded  mirror-like  ionos-  ' 

phere),  to  a difference  in  propagation  distance.  As  discussed  in  Lewis  et  al,  ^ 

the  group  delay  can  be  defined  as  the  rate  of  change  of  phase  with  frequency.  For 

the  GMH  data  presented  in  this  paper,  a finite  frequency  difference  of  1.  0 kHz  was 

used,  and  the  corresponding  phase  difference  as  a function  of  frequency  for  the 

groundwave  and  both  skywave  signals  was  obtained  by  Fourier  analysis  of  the 

respective  pulses.  The  GMH  calculations  took  into  account  ground  conductivity 

(10  mho/meter  is  assumed),  and  the  corrections  of  Wait  and  Howe  were  applied. 

Group  mirror  heights  are  plotted  as  a function  of  frequency  in  parts  C and  D 
of  Figures  3 through  11,  as  obtained  from  the  parallel  and  perpendicular  waveforms, 
respectively.  The  GMH's  are  also  presented  as  a function  of  time-of-day  for  the 
average  frequency  of  16.  5 kHz  in  figure  parts  E and  I.  As  the  frequency  of  maxi- 
mum spectral  amplitude  was  lowered  for  the  period  covered  by  this  report,  addi- 
tional GMH's  were  also  calculated  at  10.  5 kHz.  These  are  shown  in  parts  R and  V. 

The  GMH  plots  also  include  an  average  reflection  height  for  reference  purposes. 

Each  point  of  the  reference  height  is  a weekly  average,  by  time  block,  for  the 
7 -day  period  indicated.  Parts  G and  T give  t.he  average,  by  time  block,  for  the 
daily  parallel  GMH  data  of  parts  E and  R,  and  parts  K and  X give  the  corresponding 
perpendicular  GMH  averages  from  the  daily  data  of  parts  I and  V. 

4.  REFLECTION  COEFFICIENTS 

Assuming  that  the  ionosphere  acts  as  a "mirror”  at  the  GMH,  plane  wave 

7 

reflection  coefficients  were  obtained  by  comparing  the  ratio  of  the  skywave 
Fourier  amplitude  at  a specific  frequency  to  that  of  the  groundwave,  taking  into 
account  wave  spreading,  earth  curvature,  ground  conductivity,  path  lengths,  and 
antenna  patterns  Including  ground  image  effects. 

The  reflection  coefficient  ||R||  was  obtained  from  analysis  of  the  parallel  sky- 
wave  component  and  is  plotted  as  a function  of  frequency  in  part  C of  Figures  3 
through  11.  The  ||R||  coefficient  for  16  kHz  is  plotted  as  a function  of  time-of-day 
in  part  F and  ||R||for  10  kHz  appears  in  part  S.  For  comparison,  the  average  • 

6.  Wait,  J.  R. , and  Howe,  H.  H.  (1956)  Amplitude  and  Phase  Curves  for  Ground- 

wave  Propagation  in  the  Band  200  Cycles  per  Second  to  500  Kilocycles, 

Mat.  Bur.  stand.  U.'SI-Clr'c.,  p.  5f4" ^ 

7.  Budden.  K.G.  (1961)  Radio  Waves  in  the  Ionosphere,  p.  85,  Cambridge 

University  Press,  London. 
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coefficients  for  the  indicated  reference  weeks  are  also  plotted.  From  the  perpen- 
dicular skywave  pulse,  the  coefficient  n R was  obtained  and  appears  as  a function 
of  frequency  in  part  D.  The  16  kHz  ||Rj^  is  shown  along  with  its  reference  in  part  J, 
while  ||Rj_for  10  kHz  appears  in  part  W.  Parts  H and  L present  the  average,  by 
time  block,  of  the  daily  16  kHz  jjR||  and  jjRj^  data  presented  in  parts  F and  J,  re- 
spectively, while  parts  U and  Y give  the  average  of  the  10  kHz  ||R||  and  ||Rj^  data 
from  parts  S and  W. 

For  certain  coefficient  data  points,  plotted  as  asterisks(^),  the  reflection 
coefficient  appears  without  a corresponding  GMH.  For  these  particular  data,  only 
the  skywave-groundwave  ratios  could  be  obtained,  as  the  skywaves  were  too  weak 
to  provide  reliable  group  delay  information.  The  reflection  coefficients  were  there- 
fore estimated  using  a nominal  GMH  of  80  km  in  the  calculations.  These  estimated 
coefficient  values  are  included  in  the  averages  presented  in  parts  H,  L,  V and  Y, 
but  the  assumed  heights  are  not  used  in  the  GMH  averages  shown  in  parts  G,  K,  T 
and  X. 


5.  SUPPLEMENTARY  INFORMATION 

For  purposes  of  comparison  and  interpretation,  certain  supplementary  data 
are  presented.  Figure  parts  M and  N give  the  received  VLF  phase  and  amplitude 
from  the  17.  8 kHz  station  NAA  (transmitter  location:  Cutler,  Maine),  as  observed 
at  Thule  AB  over  a 3500-km  propagation  path.  Part  O of  the  figures  shows  the 
magnitude  of  the  horizontal  component  of  the  polar  magnetic  field  observed  with  a 
three-axis  fluxgate  magnetometer,  and  part  P presents  30-MHz  riometer  data,  an 
indicator  of  D-region  particle  precipitation.  These  supplementary  data  were  re- 
corded at  30-sec  intervals  by  AFGL's  Geopole  Obseirvatory  at  Thule  AB;  the  curves 
represent  the  average  of  10-min  periods.  The  solar  zenith  angle  is  given  in  part  Q 
of  Figures  3 through  11  for  the  indicated  mid-week  date. 


6.  ADDITIONAL  COMMENTS 

It  is  noted  that  a minor  Polar  Cap  Absorption  (PCA)  event  occurred  on  DAY  233. 
The  1 dB  of  signal  decrease  on  the  30  MHz  riometer  plot  of  Figure  7p  indicates 
increased  ionospheric  absorption  due  to  enhanced  energetic  particle  levels  during 
the  event.  Although  the  ionosounder  was  off  on  DAY  233,  it  was  back  in  operation 
on  DAY  234  and  the  effects  associated  with  the  event  can  be  seen  to  continue 
through  DAY  235. 
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This  paper  is  one  of  a series.  ' ' Comments  and  suggestions  for  improving 
its  usefulness  should  be  addressed  to  the  VLF/ULF  Techniques  Branch  (ETEE), 
Electromagnetic  Sciences  Division,  Deputy  for  Electronic  Technology  (RADC/ 
ETEE),  Hanscom  AFB,  Massachusetts  01731. 
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Ftfur*  3.  VLF/LF  Raflcctlvltj'  Data  for  the  Polar  lonoaphere, 
DAY  301  (30  JuIt)  - DAY  307  (36  July)  197S 
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Figur*  S.  VLF/LF  RaflvetlTlty  D»ta  for  th*  Poliir  Iono«plMr«i 
DAY  ai9  (3  Aug)  - DAY  231  (9  Aug)  1975 


Flgnr*  >.  VLF/LF  Ranectlvlty  DaU  for  the  Polar  lonosphero. 
DAY  322  <10  Aug)  > DAY  338  (18  Aug)  1978 
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Figure  7.  VLP/LF  Renectlvlty  Data  for  the  Polar  Ionosphere, 
DAY  229  (17  Aug)  - DAY  235  (23  Aug)  1975 
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Figure  9.  VLF/LF  Renecttvlty  Data  for  the  Polar  lonoephere, 
DAY  243  (31  Aug)  - DAY  249  (6  Sep)  1975 
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Flfur*  10.  VLF/LF  Rcflactlvity  Data  for  the  Polar  lonoaphere, 
MY  350  (7  Sept)  - DAY  256  (13  Sep)  1975 
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Flrare  U.  VLF/LF  R*necttvlty  »U  for  the  Polar  lonoephere, 
DAY  2S7  (14  Sept)  - DAY  263  (20  Sep)  197S 
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